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Abstract 

In this contribution we consider stochastic growth models in the 
Kardar-Parisi-Zhang universality class in 1 + 1 dimension. We discuss 
the large time distribution and processes and their dependence on the 
class on initial condition. This means that the scaling exponents do 
not uniquely determine the large time surface statistics, but one has 
to further divide into subclasses. 

Some of the fluctuation laws were first discovered in random ma- 
trix models. Moreover, the limit process for curved limit shape turned 
out to show up in a dynamical version of hermitian random matrices, 
but this analogy does not extend to the case of symmetric matrices. 
Therefore the connections between growth models and random matri- 
ces is only partial. 

1 Introduction 

In this paper we discuss results in the Kardar-Parisi-Zhang (KPZ) univer- 
sality class of stochastic growth models. We focus on the connections with 
random matrices occurring in the one-dimensional case. Consider a surface 
described by a height function x h-> h(x, t) with x G M. d denoting space and 
t 6 K being the time variable and subjected to a random dynamics. If the 
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growth mechanism is local and there is a smoothing mechanism providing 
a deterministic macroscopic growth, then the macroscopic evolution of the 
interface will be governed by 

f = v(Vh) (1) 

where u i->- v(u) is the macroscopic growth velocity as a function of the surface 
slope u. In this context, we can also focus on a mesoscopic scale where the 
random nature of the dynamics is still visible. In the famous paper of Kardar 
Parisi and Zhang [32], the smoothing mechanism is related with the surface 
tension and it takes the form uAh, while the local random dynamics enters 
as a space- time white noise rj. Moreover, the Taylor expansion of v for small 
sloped results into the KPZ equation^ 

= uAh(x, t) + \\{Vh{x, t)f + r,(x, t), (2) 

where A = v"(0) ^ in the non-linear term is responsible for lateral spread 
of the surface and lack of time reversibilitjU. 

From now on we consider the one- dimensional case, d = 1. Denote by 
h ma the limit shape, 

h ma (0 := Km ^P-. (3) 

The fluctuation exponent is 1/3, the spatial correlation exponent is of order 
2/3 [291152] . This means that the height fluctuations grow in time as 

t l/3 

and spatial correlations are t 2 ^ 3 , i.e., the rescaled height function at time t 
around the macroscopic position £ (at which h mSb is smooth)Q 

^csc (m) = m+v^t) -th^+ut-^) (4) 



1 The order and 1 in the Taylor expansion can be set to be zero by a simple change 
of (moving) frame. 

2 In more than one dimension, (V/i) 2 should be replaced by (Vft, CVh) with C a matrix. 
Then one distinguish the isotropic class, if all the eigenvalues of C have the same sign, 
and anisotropic class(es) otherwise. For instance, in d = 2 the surface fluctuation are 
very different: for anisotropic they are normal distributed in the y/ln(t) scale [2H1IS3] and 
correlation are the ones of the massless free field [51[TTj ; for isotropic is it numerical known 
that growth as t a for some a ~ 0.240 [4"8] . 

3 Whcn A = we are in the Edwards- Wilkinson class [3D] and the fluctuations are 
Gaussian with fluctuation exponent 1/4. 

4 Howevcr, depending on the initial conditions, ([2]) can produce spikes in the macro- 
scopic shape. If one looks at the surface gradient u = Vh, the spikes of h corresponds 
to shocks in u, and it is known that the shock position fluctuates on a scale t 1 / 2 . For 
particular models, properties of the shocks have been analyzed, but mostly for stationary 
growth (see [T51[TH1[H] and references therein). 
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converges to a non-trivial limit process in the t — > oo limit. Concerning cor- 
relations in space-time, it is known that along special directions the decor- 
relation occurs only on the macroscopic scale (i.e., with scaling exponent 1), 
while along any other direction the correlations are asymptotically like the 
spatial correlations at fixed time pTjlM] . The special directions are the char- 
acteristic solutions of the PDE of for the macroscopic height gradient. More 
precisely, denote by £ and r the macroscopic variables for space and time. 
Also, let 

h(£, t) : = lim r x h(£t, rt) and u(f , r) := dh(£, r)/d£. (5) 

f — >oo 

Then, u satisfies the PDE 

du/dr + a{u)du/d^ = where a(u) = —dv{u)/du (6) 

with v the macroscopic speed of growth^. The characteristic solutions of (jSJ) 
are the trajectories satisfying d^/dr = a(u) and du/dr = (see e.g. [21] 153] 
for more insights on characteristic solutions). 

The question is therefore to determine the limit process 

u i— >■ lim hr c (u) = ? (7) 

t— >oo 

One might be tempted to think that the scaling exponents are enough to 
distinguish between classes of models and therefore that the result of our 
question is independent of the initial condition. However, as we will see, this 
is not truest 

To have an intuition about the relevance of the initial condition, consider 
the fluctuations of h(0,t) with (a) deterministic initial condition, h(x, 0) = 
for all i6K, and (b) random but still macroscopically flat initial condition, 
h(x, 0) a two-sided Brownian motion with h(0, 0) = 0. For this case, the 
height function h(0,t) is correlated with a neighborhood of x = of order 
t 2//3 . In these region (at time t = 0) for (a) fluctuations are absent while in 
(b) the fluctuations on the initial condition are of order t 1//3 : this is the same 
scale as the fluctuations of h(0,t) and therefore the fluctuation laws for (a) 
and (b) will be different. 

How should we proceed to answer to our question? Literally taken is 
the KPZ equation (j2J) ill-defined, because locally one will see a Brownian 

5 In the asymmetric exclusion process explained below, one usually considers the particle 
density p instead of u. This is however just a rotation of the frame, since they are simply 
related by u = 1 — 2p. The PDE for p is the well known Burgers equation. 

6 For other observables the relevance of the initial condition was observed already in [5Tj . 
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Figure 1: Dynamics of the TASEP. Particles jump with rate one on the right 
site, but under the constraint that the site is empty. 



motion and the problem comes from the square of a white noise (in the non- 
linear term). However it is possible to give a sense of a solution of the KPZ 
equation as shown in [21131], and this solution agrees with the one coming 
from discrete approximations/models (weakly asymmetric simple exclusion 
process [7]). These works also provide an explicit solution of the finite time 
one-point distribution for an important initial condition, see jl3] for more 
explanations. 

Another point of view is to see the KPZ equation as one of the models in 
the KPZ universality class of growth models. By universality it is expected 
that the limit processes do not depend on the model in the class (but they 
depend on the type of initial condition). From this perspective, we can take 
any of the models in the KPZ class and try to obtain the large time limit. 

In the rest of the paper we consider one of such model, the totally asym- 
metric simple exclusion process (TASEP) in which the asymptotic processes 
have been unraveled. Another model for which analogues results have been 
determined is the polynuclear growth (PNG) modef]. In particular, we will 
discuss which limit distributions/processes also appears in random matrix 
theory and when the connection is only partial. 

2 TASEP 

The totally asymmetric simple exclusion process (TASEP) in continuous 
time is the simplest non-reversible interacting stochastic particle system. In 
TASEP particles are on the lattice of integers, Z, with at most one particle 
at each site (exclusion principle). The dynamics is defined as follows. Par- 
ticles jump to the neighboring right site with rate 1 provided that the site 
is empty. Jumps are independent of each other and take place after an ex- 
ponential waiting time with mean 1, which is counted from the time instant 
when the right neighbor site is empty, see Figure [TJ 



7 At least half of the limit results described below were first obtained for the PNG 
model. 
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Figure 2: Height configuration (solid line) associated with a particle (circles) 
configuration. There are three particles which can jump and the correspond- 
ing evolution of the height function is illustrated by the dashed profiles. 



Denote by t] x (t) the occupation variable of site x G Z at time t > 0, i.e., 
r) x (t) is 1 if there is a particle and if the site is empty. TASEP configu- 
rations are in bijection with the surface profile defined by setting the origin 
h(0, 0) = and the discrete height gradient to be 1 — 2r) x (t). If we denote by 
N t the number of particles which have crossed the bond to 1 during the 
time span [0, i], then the height function is given by 

f 2iV t + £* =1 (l-2 % (t)), forx>l, 
h(x,t) = I 2N t , for x = 0, (8) 

{ m - ESU+i(i - 2%(*)), for x < -1. 

as illustrated in Figure [2J 

Let us verify that TASEP belongs to the KPZ universality class. Un- 
der hydrodynamical scaling, the particle density p evolves according to the 
Burgers equation d t p + d x (p(l — p)) = 0. Thus, we have a deterministic 
limit shape. The second requirement, the locality of the growth dynamics is 
obviously satisfied. Finally, the speed of growth v of the interface is twice 
the current density, which is given by p(l — p). Being the gradient u = 1 — 2p 
it follows that v(u) = (1 — u 2 )/2, which implies v"(u) = —1^0. 

Now we discuss some of the large time results for the TASEP height 
function!!. We consider two non-random initial conditions generating a curved 
and a flat macroscopic shape. The limit processes will be called the Airy 2 
and Airy! processes, which are defined in Appendix S 

8 Most of the results has been first computed for particle positions and the statements 
described below are obtained by a geometric transformation. 
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Figure 3: Height configuration for step initial conditions. 



2.1 TASEP with step initial condition 

Consider the initial condition ^(O) = 1 for x < and r] x (0) = for x > 1, 
see Figure [3J This is called step-initial condition. 

The macroscopic limit shape for this initial condition is a parabola con- 
tinued by two straight lines: 

kmniO = < ( 9 ) 



lei, for lei > i. 



From this we have the scalin; 
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h(2u(t/2) 2 / 3 ,t) - (t/2 + u 2 (t/2) 1 / 3 ) 
hT s » := 1 ~ ' _ (t/ 2)i/3 ( 10 ) 

The large time results for the rescaled height function h r t esc are the fol- 
lowing. First, for the one-point distribution [3Tj 

lim P(^ csc (0) < s) = F 2 {s), (11) 

t— >oo 

where F2 is known as the GUE Tracy- Widom distribution, first discovered 
in random matrices [19] (see Section [3]). Moreover, concerning the joint 
distributions, it is proven [10J[TI1|32] that (in the sense of finite-dimensional 
distribution 10 !) 

lim hr c (u) = A 2 (u), (12) 

t— >oo 



9 With respect to ^ we adjusted the coefficients to avoid having them in the asymptotic 
process. 

10 In [32] Johansson the process was studied in a slightly different cut, but because of 
slow-decorrelation |17ll2~i] the present result can be proven from it. Remark also that the 
convergence in |32j is in a stronger sense. 
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Figure 4: Height configuration for a deterministic and flat initial conditions. 



where A2 is called the Airy2 process, first discovered in the PNG model by 
Prahofer and Spohn [41] (see Appendix lAll"! In particular, the Airy2 pro- 
cess is stationary, locally looks like a Brownian motion and has correlations 
decaying slow: like u~ 2 (see Figure ED- 

By universality it is expected that the Airy 2 process describes the large 
time surface statistics for initial conditions^ generating a smooth curved 
macroscopic shape for models in the KPZ class. This happens when the 
characteristic lines for space-time points on the curved limit shape go all 
together at a single point at time t — (for the TASEP and PNG are 
straight lines back to the origin). 



2.2 TASEP with flat initial condition 

The second type of non-random initial condition we discuss here is called 
flat-initial condition. In terms of TASEP particles, it is given by n x (0) = 1 
for x even and r] x (0) = for x odd, see Figure HI 

The macroscopic limit shape is very simple, /i m a(0 = \, so that the 
rescaled height function becomes 



h™(u) := "^ • (13) 



h{2ut 2 '\i) -t/2 

.£1/3 

In the large time limit, the one-point distribution of h r t esc is given byr^l 

lim P (/4 esc (0) <s) = Fx (2s), (14) 



t— >oo 



where F\ is known as the GOE Tracy- Widom distribution, first discovered 
in random matrices [5D] (see Section [3]). Moreover, as a process, it was 



11 The height function of TASEP in discrete time with parallel update and step initial 
condition is the same as the arctic line in the Aztec diamond for which the Airy2 process 
was obtained by Johansson in |33j . Extensions to process on space- like path for the PNG 
model was made in [16] . Tagged particle problem was studied in [30] , extension to space- 
like paths in [T0l|T3] and to any space-time paths except characteristic line in [17ll24| . 

12 Also for random initial conditions as in the case of Bernoulli-/}- on Z_ and Bernoulli- 
p + on Z + , p_ > p + , see [6"Pu]. 

13 For the geometric case corresponding to discrete time TASEP this result was proven 
by Baik and Rains in [5]. 
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Figure 5: Covariance Q2{u) = Cov(^4 2 (m), ^4. 2 (0)) of the Airy 2 process (dashed 
line) and gi(u) = Cov(^li('u), ^4i(0)) of the Airyi process (solid line). One 
clearly sees the difference of behavior: g 2 (w) — 2u~ 2 for u ^> 1, while g\(u) 
goes to zero super-exponentially fast. 



discovered by Sasamoto [I21II2] and it is proven that (in the sense of finite- 
dimensional distribution) 

lim h\ csc (u) = Ai(u), (15) 

t— >oo 

where A\ is called the Airy! process (see Appendix |AJ. In particular, the 
Airy! process is stationary, it behaves locally like a Brownian motion, but 
unlikely for the Airy 2 process, the decorrelations decay superexponentially 
fast (see Figure [5]), see the review [25] for more information and references. 

The Airyi process is expected to describe the large time surface behavior 
for non-random initial conditions generating a straight limit shape for models 
in the KPZ class. Unlike for the curved limit shape, the characteristic linea^l 
for space-time points for flat limit shape do not join at initial time. This fact 
is at the origin of (a) the different fluctuation behavior between curved and 
flat and (b) the difference between random flat and non-random flat. 

2.3 TASEP with stationary initial condition 

The only translation invariant stationary measure for continuous time 
TASEP are Bernoulli product measures with parameter p, p e [0, 1], which is 

14 For density 1/2 the characteristic lines are all lines parallel to the time axis. For 
density p, they have the form x = xq + (1 — 2p)t. 
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the density of particles [37J. The cases p G {0, 1} are degenerate and nothing 
happens, so consider a fixed p G (0, 1). The height function at time is a two- 
sided random walk with h(0, 0) = and P(h(x + 1,0) — h(x, 0) = 1) = 1 — p 
and P(h(x + 1, 0) - h{x, 0) = -1) = p. 

Unlike the deterministic initial conditions we need to consider the regions 
where the height function is non-trivially correlated with h(0, 0). The reason 
is that for the regions at time t which are correlated with h(at, 0), a^O, the 
dynamical fluctuations (of order t 1//3 ) and are dominated by the fluctuations 
in the initial condition (of order t 1//2 ). The correlations in TASEP are carried 
by second-class particles, which move with speed 1 — 2p. As predicted by the 
KPZ scaling, the limit 

lim M(l-2p)t + ^ 2 / 3 ,t)-(l-2p(l-p))t 

exists. The one-point distribution was derived in the PNG model with exter- 
nal source [I] and for TASEP in [27] (correctly conjectured using universality 
in [ID])- The extension to joint distributions is worked out in [3J. So-far no 
connections with random matrices for this initial condition is known. There- 
fore we do not enter in further details. 



3 Random Matrices 

We said that F\ and F2 appeared first in random matrices. We will explain 
it below and discuss whether the Airy processes also show up for random 
matrix models. 

3.1 Hermitian matrices 

The Gaussian Unitary Ensemble (GUE) of random matrices consists of Her- 
mitian matrices H of size N x N distributed according to the probability 
measured 

p GVE (H)dH = -L exp f-± Tr(if 2 )^ dH, (17) 

15 The scaling of this paper is such that the asymptotic density of eigenvalues remains 
bounded, the macroscopic variable is N and scaling exponents are easily compared with 
KPZ. In the literature there are other two standard normalization constants, which are 
just a rescaling of eigenvalues. The first one consists of replacing 1/(2N) by N: this 
is appropriate if one study spectral properties, since in the large N limit the spectrum 
remains bounded. The second is to replace 1/(2N) by 1 so that the measure does not 
depend on N: this is most appropriate if one looks at eigenvalues' minors. 
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Figure 6: Dashed line: the densities of the GUE Tracy- Widom distribution; 
solid line: the GOE Tracy- Widom distribution density. 

where dH = Y[iLidH ii Y[ 1<i< j <N dRe(H i j)dIm(H it j) is the reference mea- 
sure (and Zm the normalization constant). 

Denote by A^f ax the largest eigenvalue of a. N x N GUE matrix. Then 
Tracy and Widom in [49J proved that the asymptotic distribution of the 
(properly rescaled) largest eigenvalue is F 2 (see Figure E]): 



The parallel between GUE and TASEP with step initial condition goes 
even further. In 1962 Dyson [19] introduced a matrix-valued Ornstein- 
Uhlenbeck process which is now called Dyson's Brownian Motion (DBM). 
For hermitian matrices, GUE DBM is the stationary process on matrices 
H(t) whose evolution is governed by 



where dB(t) is a (hermitian) matrix- valued Brownian motion. More precisely, 
the entries £ M (t), 1 < % < N, Re(B id )(t) and lm(B i:j )(t), 1 < i < j < N, 
perform independent Brownian motions with variance t for diagonal terms 
and t/2 for the remaining entries. Denote by A^f ax (t) the largest eigenvalue 
at time t (when started from the stationary measure (fTTl) ). Its evolution is, 
in the large N limit, governed by the Airy 2 process: 




(18) 



dH(t) = -±-H(t)dt + dB(t) 



(19) 



lim 



A£ U m E ax (2^/3) _ 2N 



A 2 (u). 



(20) 



ATl/3 
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Thus we have seen that the connection between GUE and TASEP extends 
to the procesJ^l. 

3.2 Symmetric matrices 

The Gaussian Orthogonal Ensemble (GOE) of random matrices consists of 
symmetric matrices H of size N x N distributed according to 

p GOE (H)dH = exp Tl(H 2^ dR ^ (21) 

where dH = rii<i<j<Ar d-^ij i s the reference measure (and the normal- 
ization constant). 

Denote by A^f ax the largest eigenvalue of a N x N GOE matrix. The 
asymptotic distribution of the (properly rescaled) largest eigenvalue is F 1 
(see Figure [6]) [50] : 

(\GOE _o/\r \ 
- Jfi *b— *«)=*(«>■ (*> 

DBM is defined also for symmetric matrices: GOE DBM is the stationary 
process on matrices H(t) whose evolution is governed by 

dH(t) = ~H{t)dt + dB{t) (23) 

where dB(t) is a symmetric matrix- valued Brownian motion. More precisely, 
the entries Bij(t), 1 < % < j < N, perform independent Brownian motions 
with variance t for diagonal terms and t/2 for the remaining entries. 

We consider now the evolution of the largest eigenvalue A^ ( ^! lx (t) (when 
started from the stationary distribution f[2~Tl) ). By analogy with the GUE 
case, one might guess that in the large iV limit the limit process of a properly 
rescaled A^f ax (£) is the Airy! process. However, as shown numerically in [8], 
this is not the case. To see this, we considered the scaling (12^1) where the 
coefficients are chosen such that the variance at u = is the same as the 
variance of the Airyi process, and the covariance for \u\ <S 1 coincide at first 
order with the covariance of A\. The large time limit of the rescaled largest 
eigenvalue is denoted by B\ . Comparing the covariances as shown in Figure [7] 
we conclude that A\ ^ B\. 

hm — \ Nl/3 } =: B l{ u) * M u). (24) 
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Figure 7: Log-log plot of the rescaled correlation functions for GOE and 
GUE. For k = 1,2, we denote g k = Gov(A k {u), A k {0)) and /^ UE (w) (resp. 
f^ OE (u)) is the covariance of the GUE (resp. GOE) largest eigenvalues of a 
N x N matrix rescaled as in ( KUi) (resp. (12"1|) ). 



Thus the connection between GOE and TASEP does not extend to multi- 
time distribution^]. 



4 Conclusion 

We saw that large time fluctuations in KPZ growth models depends on the 
initial conditions. By analyzing special models we determine the limiting 
distributions and processes, which by universality should be the same for 
whole KPZ universality class. The conjecture is that when the limit 
shape is curved one gets the Airy2 process, when the limit shape is flat one 
has to further distinguish according to roughness exponent a of the initial 
condition (\h(x, 0) — h(0, 0)| ~ \ x \ a )- For a = one expects the Airyi 
process, for a = 1/2 the result from stationary initial condition. Finally, 
for < a < 1/2, at the characteristics it should still be Airy! but there 
should be a region (away of order t 1 ^ 3 ") from the characteristic lines) with a 
different (yet unknown) process, which depends on the exponent a. 

16 This connection extends partially to the evolution of minors, see |TJ[2S]. 
17 The connection extends to the top eigenvalues too [23], but not still at a fixed time 
only. 
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KPZ class Random matrices 

Curved shape Hermitian 



TASEP, step IC GUE matrices TASEP, flat IC 



t —¥ oo 



N —> oo 



F2 distribution 



A2 process 

(a) 



KPZ class 
I 

Flat shape 



t — >■ 00 



Random matrices 
I 

Symmetric 
I 

GOE matrices 
AT -> 00 



-Fi distribution 



process 



7* 

(b) 



Unknown 



Figure 8: Comparison: non-random initial conditions for KPZ vs Random 
matrices: (a) curved limit shape vs Hermitian matrices, (b) flat limit shape 
vs symmetric matrices. 



The two scheme in Figure M resume the main message of this contribution: 

(a) TASEP with step initial condition is a representant for models in the 
KPZ class with curved limit shape. For this model the large time limit 
distribution is F 2 and the limit process the Airy 2 process, A%. The 
GUE (with DBM dynamics) is a representant of random matrices with 
hermitian symmetries and also for this model the F2 distribution and 
the process A2 arises in the limit of large matrices (N — > 00). 

(b) TASEP with flat initial condition is one of the KPZ growth models with 
deterministic initial conditions and straight limit shape. The long time 
fluctuations are governed by the F\ distributions and the limit process 
is the Airyx process. On the random matrix side, GOE (with DBM 
dynamics) is a representant for symmetric random matrices and F\ is 
the distribution of its rescaled largest eigenvalue as N — > 00. However, 
it is not the Airyi process which describes its dynamical extension. 

The interested reader is referred to [28} 135] for more insights about ran- 
dom matrices and stochastic growth models. Therein a guide of literature is 
provided. Also, two summer schools lecture notes on the subject are avail- 
able [32115]. 

On the experimental side, recently Takeuchi built up a very nice exper- 
imental set-up in which the theoretical prediction (scaling exponents, one- 
point distributions and covariance) have been verified with very good agree- 
ment both for curved limit shape jH] and flat limit shape [16] . 
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A Definition of the limit processes 

In this short appendix we give the definitions of the Airyi and Airy2 processes. 
More information about these processes like properties and references can be 
found for example in the reviews 



Definition 1. The Airyi process A\ is the process with m-point joint distri- 
butions at u\ < u 2 < ■ ■ ■ < u m given by the Fredholm determinant 

m 

p ( n^ 1 ^*) - Sfe }) = det ( 1 " X«-^iX«)i»({«i,...,«m}xR) ( 25 ) 

k=l 

where Xs{uk,x) = t(x > S}.) and the kernel Ki is given by 

1 f is' - s) 2 \ 

Ki(u,s;w',s') = exp --— l(u < v!) 

^An(u' -u) V 4 ( M ~ u > J 



+ Ai(s + s' + (v? - uf ) exp \ (u' -u)(s + s') + -{v! - u) 

(26) 

Definition 2. The Airy 2 process A 2 is the process with m-point joint distri- 
butions at U\ < u 2 < ■ ■ ■ < u m given by the Fredholm determinant 

m 

p( f]{A 2 (u k ) < = det(l - XsK 2 Xs)L^{u 1 ,...,u m }xR) (27) 

k=l 

where Xs{uk,x) = t(x > Sk) and K 2 is the extended Airy kernel given by 

L dAe("'-") A Ai(x + X)M(y + A), u > u', 
K 2 (u,s;u',s') = { jR+ r (l u " (28) 

^ -/ dAe( M - u ) A Ai(x + A)Ai(|/ + A), u<u'. V ' 
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